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Abstract 



We introduce an elementary model for the electrostatic self-consistent po- 
tential in a two-dimensional electron gas. By considering the perpendicular 
degree of freedom arising from the electron tunneling out of the system plane, 
we predict a threshold carrier density above which this effect is relevant. The 
predicted value agrees remarkably well with the onset for the insulator to 
quasi-metallic transition recently observed in several experiments in Si02-Si 

and AlGaAs-GaAs heterojunctions. 
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Anderson transition |IJ in disordered solids still raises great interest among researchers. 
While metal-insulator transition (MIT) is well established in three-dimensional (3D) disor- 
dered systems, the situation is not completely understood in one- (ID) and two- dimensions 
(2D). Several decades ago, a number of papers yields the general belief that in those sys- 
tems all eigenstates are exponentially localized and a MIT no longer exists Although 
this seems to be correct in most low- dimensional systems, there exist known exceptions to 
this rule. Thus, it has been found that extended states may appear in ID random systems 
upon introducing either short-range [[SHTOll or long-range [|iT1 , [l2| correlations in the disorder. 



These purely theoretical considerations were put forward for the explanation of high con- 
ductivity of doped polyaniline || as well as transport properties of random semiconductor 
superlattices fl3fl . 

The Anderson transition induced by diagonal disorder at the band center in finite 2D 



systems was already studied by Yoshino and Okazaki |pB) . The relevance of these midgap 
states for the metallic conductance of 2D systems has been discussed in detail by Liccia- 
rdello and Thouless [|T5[]. These authors suggested that between the mobility edges there is 
a tendency for the conductance to decrease slowly as the sample size is increased and that 
they may be no absolute minimum metallic conductance. However, recent experiments have 
provided clear evidences of the MIT-like in high-quality 2D electron and hole systems. After 



the pioneering work by Kravchenko et al. on electrons in Si ||T6Hlq] and more recently by 



Hanein et al. on holes in GaAs [I5H, it has be become clear that 2D gases undergo a crossover 



from an insulating regime at low density to a metalliclike behavior at high density, where 
the quasi- metallic phase is characterized by a strong decrease of the resistivity as the tem- 
perature decreases. Discrepancies between the standard one-parameter scaling theory 0, 
establishing that the 2D gas should be insulating, and the above mentioned experiments are 
usually attributed to a strong electron-electron interaction. A number of theoretical models 
have then been proposed to explain the observed transition in 2D systems, ranging from 
new liquid phases |2U|| , Wigner glass ||21|| , spin-orbit induced transition |22|| , decoherence 



due to quantum fluctuations f23|, superconducting phase |24| and anyon superconducting 



model ||25|| . The basic ingredient of these models is the assumption that the system is purely 



2D and, consequently, new phenomena are to be considered to explain the observed tran- 
sition. In this work we undertake a different way by considering the electronic motion in 
the perpendicular direction. In so doing, we find the conditions under which this degree of 
freedom could be relevant. Surprisingly, the critical carrier density leading to perpendicular 
motion is close to that determined in experiments to observe the transition. Our main aim 
is then to point out the importance of this perpendicular degree of freedom, which should 
be taken into account in more elaborated models. 

Our approach is based on the competition between the confining potential, appearing 
at the hetero junction even at zero gate voltage, and the repulsive potential arising from 
the excess carrier at nonzero gate voltage. To proceed, let us start by considering the 
hetero junction at zero gate voltage. Since we are only interested in the basic phenomena 
without entering in many details, we make use of a simple variational Hartree calculation 



presented in standard textbooks ||26|| . This will make our reasonings clearer while keeping 
a good qualitative description of the involved physics. The envelope function of the lowest 
subband of a 2D electron gas is reasonably well accounted for by the Fang-Howard trial 
function: 

X ^ = {^2 Z exp ir\ bz ) ' ^ 

where z is the coordinate along the growth direction (perpendicular to the heteroj unction) 
and b is determined by minimizing the kinetic energy plus the Hartree energy per electron 



(see Ref. [26[ for details). In terms of the effective Bohr radius, a*, and the 2D electron 
density, n°, the variational parameter is roughly given by 

/16vm°\ 1/3 



V a* 



(2) 



while the maximum value of the Hartree potential is |26 



e 2 n° 



^o = 3-r* (3) 
eb 
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e being the dielectric constant of the medium. The Hartree potential increases smoothly 
from zero up to Vq on increasing the distance z > from the heteroj unction. Thus, at zero 
gate voltage the electrons lie on the lower subband and confined, thus forming a 2D gas. 
The lower subband energy is expressed as: 

£l = 2.5^. (4) 
eb 

The situation may be different when the gate voltage induces an excess carrier density 
An s , as we will show below. This excess carrier is not compensated by any other charge 
close to the heterojunction, thus leading to a local negative charge density. As a crude and 
first approximation, we assume that this excess charge density is confined to a plane close 
to the heterojunction. By solving the Poisson equation, one can obtain the potential energy 

due to this charged plane, namely —eFz for z > 0, where 

_ eArig 

2e v ; 

Due to this repulsive potential, electrons may tunnel through the barrier formed by the 
Hartree potential plus —eFz and escape. To account for the tunneling process, we replace 
the actual potential by that depicted in Fig. [I]. This replacement is not essential in the 
calculations since both the confining Hartree potential and the repulsive potential due to 
the charged plane are known, but it allows us to obtain closed analytical expressions for 
the transmission coefficient. Since the Fermi energy is usually very small in this system, we 
can also safely neglect it when calculating the transmission coefficient and take E — ei, as 
indicated in Fig. [3]. Thus, the classical turning points are Z\ and z 2 = Zi+(Vo— e\)/eF, shown 
in the figure. From the WKB semi-classical approximation, the transmission coefficient is 
given by 



T = exp 



exp 



2m* f Z2 I 

v ~W J dz v v ° ~ £l ~ eF{yZ ~ Zl "' 



4 [2m* (V Q - £i) 3/2 



(6) 



3 V h 2 eF 

The onset for the crossover from a 2D to 3D behavior can be determined from the condition 
that the transmission probability is large. We then assume that this crossover appears when 



the exponent in is of the order of unity in absolute value. This occurs for a critical 
density n c s = An s + n° s such that 



4 ^w* {y,- ei f/ 2 

3V1F^F (7) 

Since F is a function of An s , we can readily determine the value of the critical density from 
the above condition for which the probability for tunneling out of the well is large. Using 
(f|) and (f|) we finally obtain n c s ~ (5/3)n°. Since the carrier density at zero gate voltage is 
related to the Fermi energy E F by the expression n° = 2m*E F /h, we finally arrive at the 
condition 

n c s ~ 7E° F m* x 10 11 cm" 2 , (8) 

where the Fermi energy is measured in meV and the effective mass in units of the free 
electron mass. 

Now let us compare our qualitative prediction @ with the experimental values. A 2D 
electron gas in Si has been studied by Kravchenko et al. ||18|| , who observed the transition 
for an electron density 0.85 x 10 11 cm~ 2 . The Fermi energy was Ep = 0.6 meV and the 
effective mass m* = 0.2m. Inserting both values in (^|) we obtain a critical density n° s ~ 
0.84 x 10 11 cm -2 . On the other side, 2D hole gas in GaAs was demonstrated to undergo a MIT 



at a hole density 0.15 x 10 11 cm~ 2 by Hanein et al. [19|. Taking the values Ep = 0.04 meV 
and m* = 0.4m we get from (|j) that n c s ~ 0.11 x 10 11 cm -2 . The agreement should be 
regarded as surprisingly good, in view of the crude approximation we made to obtain it. 
As a conclusion, our model points out the relevance, under some circumstances, of the 
perpendicular degrees of freedom in the so-called 2D electron gases. As soon as the electron 
gas becomes a non-perfect 2D system, the scaling theories predicts the occurrence of a MIT 
transition like that recently observed. 
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FIGURES 




Slope=-eF 



FIG. 1. Schematic diagram of the band-edge profile of the heterojunction at high carrier density. 
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